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Abstract. We develop an analytical scattering formalism for computing the 
transmittance through periodic defect lines within the tight-binding model of graphene. 
We first illustrate the method with a relatively simple case, the pentagon only defect 
line. Afterwards, more complex defect lines are treated, namely the zz(558) and the 
zz{5757) ones. The formalism developed, only uses simple tight-binding concepts, 
reducing the problem to matrix manipulations which can be easily worked out by any 
computational algebraic calculator. 

PACS numbers: 81.05.ue, 72.80.Vp 
1. Introduction 

Grain boundaries (GBs) in artificially grown solids are, most likely, unavoidable. This 
is particularly true for solids grown by chemical vapor deposition. In this method, 
the crystal starts growing simultaneously at different locations on the substrate. The 
relative orientations of the domains have a stochastic distribution and when two domains 
growing at different locations approach each other they form a GB. 

If the grown crystal is used in a device larger than the size of the grains, then an 
electron has to pass through one or several GBs as it makes its way across the device. 
Therefore, the scattering problem of an electron off a GB becomes technologically 
relevant. 

In low dimensional systems, such as graphene [1], the study of GBs is an active 
field of research [21 E] and it has been shown that this type of disorder has a strong 
impact in the transport properties of graphene [UEIE]. Indeed, one-dimensional defect 
lines can give rise to exotic effects, such as a valley filter [71 [8], where states from one 
Dirac cone are filtered from those belonging to the other inequivalent cone. It has 
been suggested that GBs can be exploited for novel graphene-based nanomaterials and 
functional devices [6J. 
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The study of scattering by extended defects is revealing itself of increasing interest, 
specially after the recent work [9] by Tsen et ai. In the latter, using chemical vapor 
deposited poly crystalline graphene, the authors made electric measurements across a 
single GB. In this way, they were able to measure the electronic properties of single 
GBs. Tsen et al. have found that the transport properties of these systems are strongly 
dependent on their microscopic details. Each CVD synthesis method, typically gives rise 
to GBs with similar resistivity profiles. But the GBs originating from different synthesis 
procedures are normally very distinct. All this highlights the importance of the studies 
presented on this paper concerning the electronic properties of defect lines in graphene. 

The study of one localized impurity in a one-dimensional system has a simple 
analytical solution [10] , which was later generalized to the ladder case [11] . This method 
is, however, inappropriate for tackling more complex systems. For problems of the same 
nature as the ones considered in this work, the method of non-equilibrium Green's 
functions is often employed [TJ [12] . Unfortunately, such method is not so easy to follow 
by the non-specialist, although there is a one-to-one correspondence between the Green's 
function method and the mode matching one [13]. A more elementary method, based on 
the work of Ando [H] , was applied to the solution of a defect on ultra-narrow graphene 
|15j . Here we develop an approach, with close resemblance to the method of Ando for 
quantum point contacts, which is particularly suitable to deal with extended periodic 
defects. Its starting point is the reduction of the two-dimensional scattering problem to 
an effective quasi-one-dimensional one, whose solution is simpler to work out. 

In the context of the first-neighbor tight-binding model of graphene, we study 
the electronic scattering from periodic defect line in graphene, developing a systematic 
procedure to approach such problems. We illustrate the method with three types of 
defect lines oriented along the zigzag direction. In a companion paper we have focused 
on the continuum low-energy limit of the electronic scattering from the defect lines 
under study in the present text [16] . 

In Section [2] we start by describing briefly and in general terms how these kind 
of problems can be tackled. After that introduction, we concentrate on the study 
of electron scattering from a simple model of a defect line on graphene, namely the 
pentagon- only defect line (Section [s]). We then analyze two more complex defect lines, 
namely the zz{h7b7) [16] and the zz{bbS) ones [IZl El [71 HE] (see Section |4]). Finally, in 
Section |5] we summarize the results obtained and pinpoint the strengths of the method 
presented in this paper. 

2. General formulation of the problem 

Electronic scattering from periodic extended defects in the center of a pristine 2D crystal 
is essentially equivalent to the electron scattering from a localized defect at the center 
of a quasi-lD crystal. Such a statement arises from the fact that the periodic extended 
defect preserves crystal's translation invariance along the defect direction. Therefore, 
we can Fourier transform the 2D crystal along the latter direction thus converting the 
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2D problem into a quasi-lD equivalent one. 

Thereupon, this class of 2D problems can be treated under the framework of 
scattering ID problems. Accordingly, in the present section we describe, with some 
generality, how one can work out electron scattering from a defect located at the center 
of a quasi-lD crystal modeled by a first neighbor tight-binding (TB) model (see Fig. 
[T|. This will be the starting point for the treatment of some 2D scattering problems in 
graphene (see Sections IS] and lij) . 
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Figure 1. (Color online) Scheme of a general quasi-lD crystal with a defect at position 
n = 0. The unit cell in the bulk (green rectangles), encompasses r Wannicr states, 
while in the defect (yellow circumference) there are x atoms. The generalized hopping 
amplitude Vl {Vr — V^), is a matrix that contains all the hopping amplitudes between 
atoms of neighbor bulk unit cells. The generalized hopping amplitude Dl {Dr), is a 
X X r matrix containing the hopping amplitudes connecting all the atoms in the defect 
and those in the unit cell at n = — 1 (n=l). 



We start by writing the TB bulk equations of the general quasi-lD crystal with a 
defect at its center (see Fig. [T]). Away from the defect, these can be cast in the form 

VRc{n -1) + {H- elr)c{n) + VLc{n + 1) = 0. (1) 

In Eq. ([T]), c stands for a column vector with as many entries as there are Wannier 
states in the quasi-lD crystal's unit cell; we denote this number by r. Therefore, the 
terms H,Ir, Vl and Vr are r x r matrices, 1^ standing for the unit matrix. The matrix 
H describes the hopping processes occurring inside the unit cell, while Vl (V/j) describes 
the hopping processes occurring between the unit cell at position n and the unit cell at 
position n + l [n — l). The hermitian nature of the Hamiltonian requires that Vr = vl 
The TB bulk equation [Eqs. Q] can usually be written in a different form, where 
the amplitudes of the unit cell located at position n + l are expressed in terms of the 
amplitudes of the unit cells located at positions n and n — l, 

c(n + l) =Tic(n) + T2c(n-l), (2) 

where Tj are r x r matrices, which we call transfer matrices. These matrices are generally 
non-hermitian. Usually Eq. (|2| cannot be obtained directly from Eq. ([T]) for site n, 
because the matrix Vl is not invertible. However, for the cases of interest, we can in 
general obtain Eq. ^ from the TB equations [Eq. ([T|] for the sites n and n + l. When 
the one-dimensional chain is ^i?-like (see Fig. |2| the transfer matrix description further 
simplifies to 

c(n+ 1) = Tc(n), (3) 
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where there is now only one transfer matrix T relating neighbor amplitudes. Since 
in the present text we are specially interested in studying graphene crystals, whose 
Fourier transformed systems give rise to AB-like chains, we will from now on simplify 
our analysis by assuming that our general ID chain is described by a relation of the 
form of Eq. ^ . 
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Figure 2. (Color online) Scheme of a general quasi-lD AB crystal. The unit cell of 
the crystal (green rectangles) encompasses r = a + b Wannier states. The generalized 
hopping amplitude Vl {Vr = Vl) , is a matrix that contains all the hopping amplitudes 
between atoms of neighboring unit cells. These generalized hopping amplitudes, are 
only non-zero when connecting the a and b atoms of different unit cells. In contrast, 
they are zero both between the a atoms of two neighboring unit cells, and between b 
atoms of neighboring unit cells. 

Inherent to the study of any scattering problem is the determination of the 
transverse propagating (or evanescent) modes allowed in the system. In the case of a 
one-dimensional periodic chain the computation of the modes can be done using Bloch's 
theorem. In that prescription, we write the amplitudes at position n+l and n—1 in terms 
of those at position n: c(n+l) = e''^^"'c{n) = \jc{n) andc(r2— 1) = e"'^J"c(n) = c{n)/\j, 
where kj stands for the wave-number along the chain direction associated with the mode 
indexed by j, and a stands for the length of the primitive vector. Using these relations, 
we can rewrite Eq. ([T]) as 

^VRc{n) + {H- elr)c{n) + A,Vic(n) = 0. (4) 

The r eigenvalues \j = e*'^^" and their r associated eigenvectors permitted at a 
given energy e, are obtained by solving this generalized eigenproblem. On the other 
hand, when one is able to write the TB equations in the form of a recurrence relation 
of the form of Eq. ([s]), it becomes obvious that the \j and \^j) are, respectively, the 
eigenvalues and right eigenvectors of the transfer matrix, T. 

Let us now turn our attention to the defect, so that we can determine the boundary 
condition that the modes should satisfy when scattering from it. We start by assuming 
that regardless of the complexity of the defect at the center of the ID crystal, the latter 
is located at position n = (see Fig. [T]). Therefore, the TB equations associated with 
the amplitudes at the defect can be generally written as 

DRci-1) + (Hd - l..)d(0) + Dic(l) = 0. (5) 

In the above equation, Dl (Dr) is a. x x r rectangular matrix containing all possible 
hopping amplitudes connecting the defect Wannier states and those of the unit cell 
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located at position n = +1 [n = —1). The matrix Hd encompasses all the hopping 
amplitudes between the x Wannier states of the defect, whose amplitudes are grouped 
in the x-dimensional array d(0). Note that Eq. Q is a shorthand for a system of 
equations, containing x TB equations associated with the x states at the defect and 
X + 2r unknown amplitudes. As a consequence, if we want to write a passage equation 
relating the amplitudes at c(l) and those at c(— 1), we will also have to use the TB 
equations at position n = —1 

Vrc{-2) + {H- I.)c(-1) + DLd{0) = 0, (6) 

as well as the transfer matrix relation c(— 1) = Tc(— 2). Doing so, we end up with x-\-3r 
unknown amplitudes and x + 2r equations. We choose to solve them by expressing 
c(l), d(0) and c(— 2) {x + 2r amplitudes) in terms of c(— 1). Therefore, we will be able 
to write a passage equation relating the amplitudes before and after the defect as 

c(l) = Mc(-l), (7) 

where M is a r x r matrix. It is natural to interpret the latter matrix as a boundary 
condition matrix imposed by the defect on the wave-function at each one of its sides. 

We have so far shown how one can determine both the scattering modes of the one- 
dimensional chain and the boundary condition they must obey at the defect. The only 
step remaining is the computation of the scattering coefficients. Let us suppose that we 
have an incoming mode from n = — oo. We then expect to have both transmitted modes 
outgoing to n = +oo and reflected modes outgoing to n = — oo. As a consequence, we 
can write the wave-function at each side of the defect as 

\^in) = \^^in)) ^Ar+^|^>) + X]p,,A^+V/), forn < -1 (8a) 

\^{n)) = \^\n)) = J2^^J\];'\i^f), forn > 1 (8b) 

i=i 

where the superscript L stands for the wave-function to the left of the defect {n < 0), 
while the superscript R stands for the wave- function to the right of the defect {n > 0). 
In Eqs. ([S]), the sum over j up to r> (r<) stands for a sum over all the p> (p<) 
propagating modes moving in the direction of n = +oo {n = — oo) and all the s> 
(s<) evanescent modes, decreasing (increasing) with increasing n. Note that the total 
number of propagating and evanescent modes must add up to r, the total number of 
modes: r = r>+r< =p> + s>+p< + s<, which in the AB chain happens to be equal 
to the number of atoms in the unit cell. Moreover, the coefficients pij (Tij) stand for the 
reflection (transmission) coefficients of the modes indexed by j for an incoming mode i. 

Before proceeding, let us emphasize that in this text, we are going to use 
interchangeably two types of notation: we will use the vectorial notation when referring 
to the arrays of amplitudes at a given unit cell located at position n [see Eqs. ([T])-( 
[7|]; the Dirac ket notation will be used when referring to the eigenmodes of the one- 
dimensional chain, as well as to the wave-function at a given position n [see Eqs. 
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As long as we know the mathematical expressions for the modes of the quasi- 
one-dimensional chain, \ipf) and \ipf), and the boundary condition matrix M, the 
determination of the scattering coefficients in Eqs. ([s]), pij and r^j, is a straightforward 
calculation. 

If we define a matrix f/, with columns which are the r transverse modes 



U 



l^>),...,ie),lv^<),---,IO 
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it is clear that 

|*^(r2)) = f/|$^(n)), 
l^^(n)) = f/|$^(n)). 

The passage equation, [7], then becomes 

|$«(1)) = f/-iMf/|$^(-l)). 

This is a non-homogeneous system of r linear equations with r unknowns, pn 
and Tji, • • • , Tjj.^, which we solve to obtain these scattering amplitudes. 

When carrying out these calculations in a computer algebra system, the following 
remark may be useful. Since the transfer matrix is non-hermitian, its eigenbasis is not 
orthogonal. Nevertheless we can define a dual basis {ipj), j = 1, . . . ,r hj the relation 



(11a) 
(lib) 

(12) 

• ) Pir< 



and it is clear that the m atrix U ^ j^^ row is just the vector {ip 



(13) 

seen that the i^^ 

column of U is in the Appendix A, we show that these vectors {ipjl can be simply 
obtained as the right eigenvectors of the transpose of the transfer matrix T"^. 

It is worth noting that the boundary condition arising from the presence of the 
defect, Eq. ([t]), must conserve the particle current (see Appendix B). Equivalently, the 
current on the left hand side of the defect, JTl = i7(— 1) = ("^(—1)1 J7'|'^(— 1)), niust be 
equal to the current on the right hand side of the defect, J'r = J^{1) = ('?/'(l)| i7|^(l)), 
where stands for the current operator. Therefore, the boundary condition matrix, M, 
must satisfy the following equality 



J. 



(14) 



In the following sections we solve three electronic scattering problems using the 
machinery just presented. We will see that it simplifies the mathematical treatment of 
such problems without leading to a substantial loss of physical insight and intuition over 
the physical phenomena going on. On one hand, using the transfer matrix formalism 
we are able to write the scattering modes = e'^''^°'\ipj) directly in terms of the 

energy e and the wave-number along the defect line, kr^. We thus avoid to write the 
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former in terms of the wave- number perpendicular to the defect hne, k ■ U2. On the 
other hand, the described procedure allows us to compute an expression relating the TB 
amplitudes at the two sides of the defect. This has the obvious interpretation as being 
the boundary condition imposed on the wave- function at the defect. As we will see in 
Sections |3] and |4| such a boundary condition relation can be obtained after some simple 
algebraic manipulations of the TB equations at the defect. 

With this in mind, let us now apply this formalism to three types of zigzag oriented 
defect lines in graphene: the pentagon-only defect line (see Fig. |3]), the zz{5757) defect 
line (see left panel of Fig. [?]) and the zz{55S) defect line (see right panel of Fig. [?]). 

3. Scattering from a pentagon- only defect line: a pedagogical example 

In what follows, we shall illustrate the procedure described in Section |2] by working out 
the scattering of electrons from a pentagon- only defect line in the context of the first 
neighbor tight-binding model of graphene (see Fig. |3]). As this system is invariant under 
translations along the defect line, Bloch's theorem allows the Fourier transformation 
along that direction. Therefore, the 2D graphene layer can be transformed into an 
effective ID chain (see Fig. |4]). This effective ID chain depends on the quantum number 
kx, the wave- number along Ui direction. Based on this ID effective description, we can 
straightforwardly employ the procedure described in Section |2} 

3.1. Formulating the problem 

Let us assume we have the following system: a defect line of the type represented in 
Fig. |3} in an otherwise perfect graphene lattice. 

As we can see by inspection of Fig. [3} the lattice vectors ui and U2, in cartesian 
coordinates, read 

Ui = a(l, 0), (15a) 

U2 = a(-^,^), (15b) 

where a stands for the lattice parameter. 

Using this coordinate system, we can write the first neighbor tight-binding 
Hamiltonian of the system sketched in Fig. [s] as -ff = H'^ + + H^, where H'^ 
(H^) stands for the Hamiltonian above (below) the defect line, while the remaining 
term, H^, describes the defect line itself. In second quantization the explicit forms of 

and read 



6^(m, n — 1] 



-f 6^(m — 1, n — 1) d(m,n) + h.c.j, (16) 
where for H'^ (H^), n > 1 {n < —1). Moreover, reads 

= -t^|[^d"f(m + l) + a"^(m,0) + 6"f(m,0) d(m) + /i.c.}, (17) 
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Figure 3. (Color online) Graphene's honeycomb lattice with a zigzag oriented linear 
defect, which we dub pentagon-only defect line. In the bulk the hopping is t, while 
between atoms of the defect it is ^t. The lattice vectors are denoted by Ui — a(l,0) 
and U2 = (— l,Vl3)a/2, where a stands for the lattice parameter. 



where t is the usual hopping amplitude of pristine graphene, while ^t is the hopping 
amplitude between the atoms of the defect line, as represented in Fig. |3] 

As was referred in the beginning of this section, this system is invariant under 
translations r = mui (where m is an integer). Thus, we can make use of Bloch's 
theorem and Fourier transform the Hamiltonian along ui . This diagonalizes the problem 
relatively to the variable m, introducing a new quantum number, k^. Consequently, the 
resulting Hamiltonian turns out to be equivalent to the Hamiltonian of a one-dimensional 
chain with two atoms per unit cell and a localized defect at its center (see Fig. |4]). 

The Hamiltonian of the effective one-dimensional chain, can be written as H{kx) = 
H^{kx) + H^{kx) + H^{kx). Its three terms read 



H^'ikx) 



n 



\^'b\kx., n — 1) + tP{kx, n) a{kx, n) + /i.c.|, 
2^tcos{kxa)d^{kx)d{kx) — td\kx,0)d{kx) 



;i8a) 



+ tb\kx,0)d{kx) + h.c. . (18b) 

This one-dimensional chain has alternating hopping amplitudes between the atoms, t 
and f = t(l + e*'^"='^). At the defect, there is a on-site energy term, 'e{kx) = —2^t cos{kxa) , 
which depends on the value of the longitudinal momentum k^- 

The TB equations of such system are, as usually, obtained from Schrdinger's 
equation 

H{kx) l/i, kx) = e^_fc, l/i, kx) . (19) 
In Eq. (19) H{kx) stands for the Hamiltonian of the effective one- dimensional chain. 
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Figure 4. (Color online) Effective one-dimensional chain obtained after Fourier 
transforming the Hamiltonian of graphene with a pentagon-only defect line, along 
the Ui direction. 



while the eigenstate k^) can be expressed as a hnear combination of the site 
amphtudes along the one- dimensional chain 

+ 00 

= ^ ^A{k^,n)\a;k^,n) + B{kx,n)\b;k^,n) + D{k^)\d; k^), (20) 

n=—oo 

where the |c; k^, n) = c^lk^, n)\0) stands for the one-particle states at the atom c = a,b,d 
of unit cell n of the one-dimensional chain. 



3.2. Bulk properties 

From what we have written above, it is simple to conclude that the TB equations in the 
bulk (n 7^ 0) of the one- dimensional chain (see Fig. |4]), read 

eA{k,, n)= - tB{k,, n) - {t'yB{k„ n-1), (21a) 
eB{ka:,n)= - tA{k^,n) - t' A{k^,n + I). (21b) 

To recast these equations in the form of a transfer matrix relation, we solve Eq. 



( |21b| ) for A(A;^, n + 1) and Eq. ( pla| ) for B{k^, n 
A{k^,n + 1) = 



B{k^,n) 



j^A{k^,n) - ^B{k^,n). 



A{kr,,n) 



t' 
it') 



-B{k^,n 



(22a) 
(22b) 



where we have used t' = t(l + e """). If now we write Eq. (22b) for n + 1 and substitute 



Eq. (22a) in it, we obtain a recurrence relation between the amplitudes at unit cell n + 1 



and those at unit cell n 

L{n + l) = T{e,k^).L{n), 



(23) 



where L(n) = [A{kx^n)^B(kx^n)] . The passage matrix T{e,kx) has the explicit form 



T(e, k:r. 



2 COS ( ^ 



L f 
■f 4cos2(^ 



t2 



(24) 



As said in Section |2} the eigenvalues and eigenvectors in the bulk are obtained from 
the diagonalization of the matrix T{e,kx). In particular, eigenvalues of the transfer 
matrix which have unit modulus, |Ap = 1, correspond to Bloch solutions propagating 
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along the one-dimensional chain (a band state); eigenvalues with a modulus that is 
different from one, |Ap 7^ 1, correspond to evanescent Bloch solutions. These evanescent 
solutions decrease with n — )■ +00 {n — > —00) when |Ap < 1 (|Ap > 1). 



3.3. The defect 



We now go on to write the TB equations at the defect (see Fig. |4]), and work them out 
in such a way that we can write a boundary condition in the form given by Eq. ([T]). 
The TB equations at the defect (see Fig. |4j)are: 



eA{k,, 1) = 


-{t'rB{k^,0)-tB{k,,iy, 


(25a) 


eB{K,0) = 


-tD{k^)-t'A{k^,iy, 


(25b) 




- t{A{k,, 0) + B{k,, 0)) - 2^t cosik,a)D{k,); 


(25c) 


eA{K, 0) = 


-{tyB{k,,-l)-tD{k,); 


(25d) 


eB{k,,-l) = 


-t'A{k,,0)-tA{k,,-l). 


(25e) 



As we have mentioned above, the aim is to obtain an equation relating the wave 
function at the two sides of the defect. With that in mind, we solve each equation in 
(25) for the rightmost amplitude appearing in it: 

(tr 



B{k,A) 
A{k^A) 

Bik,,0) 

D{k,) 

A{k,,0) 



t 



A{kx, 1 J 



t 



Bikx,0), 



Dik,) - Aikx,0), 



^B{kx,0)-^-D{kx) 

e + 2^t coiiikxa) 
t 

^A(A:,,0)-^i?(A;,,-l), 
^B{kx,-l)-^-A{kx,-l). 



It is convenient to write the above set of equations in matrix form as 



B{h,,l) 
A{kx, 1) 

A{h,A) ' 
B{k,,0) 

5(fc,,0) " 
Dik,^ 

D{k^) 
A(A;,,0) 

A(A;„0) 
B{kx,-l) 



ijr 
t 



i 



t' t> 

-1 



A{kxA) 

B{k,,0) 

Bikx,0) 
D{kx) 



+ 2^ cos{kxa) 1 






-1 


e 


(t'Y 
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-1 
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t' 


t' 


-1 






Dik,) 
Aikx,0) 



Aikx,0) 
B{k,,-1) 

B{k,,-1) 
A{k,,-1) 



(26a) 
(26b) 
(26c) 
(26d) 
(26e) 

(27a) 
(27b) 
(27c) 
(27d) 
(27e) 
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It is now straightforward to write the boundary condition connecting the two sides 
of the defect as 

L(l) =M55.L(-1), (28) 

where the matrix M55 is a 2 x 2 matrix defined by 

M55 = R.Niie, h).N2{e, h).N3{e, h).Ni{e, h).N2{e, h).R^; (29) 

the matrix R is the Pauh matrix, used to interchange rows A and B. The matrices 
Ni, N2, and N3, after substituting t' = t{l + e*''""), read 

(1 + e-''"-") 



Ni(e,fc,.) 
N2(e,A;,) 
N3(e,A;,.) 



t 

-1 







]_ _|_ gifei^a 



f 1 



e+2t^ cos(kxa) 
t 

-1 







(30a) 
(30b) 
(30c) 



From Eqs. (28)- (30) we conclude that matrix M55, and consequently the boundary 



condition imposed by the defect, depend both on the energy, e, and on the longitudinal 



momentum, k^. 



3.4- Computing the scattering coefficients 

Using the spectrum of modes allowed in the bulk of the one-dimensional chain [obtained 
from the transfer matrix T(e, k^)] and the boundary condition calculated above [see Eqs. 



(28)-(30)], we can now solve completely the scattering problem from a pentagon-only 



defect line. 

We consider the scattering process in which we have an incoming mode from 
n = —00. In this situation, due to the presence of the defect at n = 0, there will 
be a reflected mode as well as a transmitted one. For now, let us suppose to be working 
at positive energy and around the = +1 Dirac point [K+ = (47r/3a, 0)]. 

We choose to denote |\E'>) (|\E'<)) as the positively (negatively) moving mode of the 
transfer matrix when we are around the K_|_ Dirac point. Furthermore, the eigenvalue 
associated with this mode is going to be denoted by A> (A<). Therefore, we can write 
the wave-function on each side of the defect as 



|^(n < 0)) = A^+^|^>) + pA^+^|^<), (31a) 

|^(n > 0)) = rA!;-i|^>), (31b) 

where p and r stand for the reflection and transmission scattering amplitudes, 
respectively. 

The direction of propagation of each of the modes obtained by diagonalization of 
the transfer matrix, can be determined by calculating the corresponding current. It is 



shown in Appendix B that, for graphene, the current operator along the equivalent ID- 



chain is tay/h; in Fig. Bl we plot the current associated with each one of the transfer 
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matrix's modes. For different clioices of e and the incoming, reflected and transmitted 
modes must be cliosen accordingly with the direction of propagation of the current. 

In order to determine the scattering coefficients p and r, following the procedure 
detailed in Section |2| we write Eq. (12) for this specific case, assuming an ordering of 
the modes as {|\t'>), |^<)}; the resulting equation is 



T 





u- 



^55 



u 



n.=+l 



1 
P 



Solving the above linear system is trivial, resulting in 



P = - 



'21 



T 



det {U-^l 



^55^^)22' 



{u-m,,u) 



(32) 

(33a) 
(33b) 



22 



3.5. The transmittance 

In the context of graphene, we are mostly interested in the low-energy region of the 
spectrum, that is close to the Dirac points. Therefore, in what concerns the electronic 
transport, it is a natural choice to plot the transmittance, T = |rp in terms of 

q 



k — where k is measured from the center of the zone [see Fig. |5]. 
In the low-energy limit, the graphene quasi-particles impinging on the defect line are 
massless Dirac fermions with a wave-vector q, which defines their propagation direction. 
Therefore, in this limit, it is intuitive to refer the transmittance to the angle 6 that 
q makes with the defect line. As a consequence, and despite the fact that we are 
not necessarily working at low energies, we will choose to express the transmittance 
(obtained from the TB model) in terms of e and 6, instead of doing so in terms of e and 



kx- In Appendix D we show how starting from the expressions of T(e, k^) we can obtain 
the transmittance in terms of the energy and the angle 6, namely T(e, 6). 



In Fig. |B1| it can be seen that the two modes of the transfer matrix have opposite 
directions of propagation. The figure also shows that when the sign of the energy is 
changed, the direction of propagation of the modes is reversed. Suppose that for a 
given pair (e, kx), the mode identified by 9 (or +qy) is the one with positive direction of 
propagation. Therefore, if we make the change e — )■ — e, then the mode propagating in 
the positive direction along the one- dimensional chain is now the one with —6 (or —qy) 
[compare panels (bl) and (b2) of Fig. [sj. 

As a consequence, and for the sake of comparison between scattering processes 
occurring at positive energy and those happening at negative energy, we will plot the 
transmittance against the angle 6' = 6 when e > 0, and against the angle 6' = —9 when 
e < 0. 

In Fig. |6]we represent the transmittance T = |rp as function of the angle 9 (when 
e > and —9 when e < 0), made between q of the incoming particle and the barrier [see 
Fig. [5] . In the different panels of Fig. |6| we present the transmittance curves for several 
values of the energy. The curves in the latter figure refer to the Dirac point K+. Those 
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Figure 5. (Color online) (a) Pristine graphene FBZ where the incident wave vector, 
k, is signaled, as well as the vector q = k — K^,. (b) Schemes of the electron scattering 
across a defect line in 2D for low energy (wave vector around the Dirac cone) in two 
different cases: (bl) e > 0; (b2) e < 0. The labels (i), (r) and (t) stand respectively 
for the incident, reflected and transmitted modes' wave-vectors. 



referring to the Dirac point K_, are mirror symmetric to the former ones relatively to 
the axis 9 = n/2. 



e/n e/7t 




e/ji e/7t 



Figure 6. (Color online) Transmittance across the pentagon- only defect line as 
function of the incoming angle 9 between q (of the incoming particle) and the barrier. 
Remember that for negative energies, the transmittance is plotted against the angle 
—6. The hopping parameter at the defect, ^, was fixed to ^ = 1 in this set of plots. 
The transmittance is plotted for several different modulus of the energy. These are: 
(a) \e/t\ = 0.01; (b) |e/i| = 0.05; (c) \e/t\ = 0.1; (d) \e/t\ = 0.5. Positive energies 
are represented by the black full lines, while the negative ones are represented by the 
dashed red lines. Only the curves associated with the Dirac point K+ are represented. 
Those for the Dirac point K_ are obtained from the former by a reflection of these 
over the axis 9 — Tr/2. 
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Modifications of tlie flopping parameter C, at the pentagon- only defect resuft in very 
different transmittance befiaviors. Tfie iow-energy resuft is in good accordance with 
that obtained in Ref. [16] using the continuum approximation. 



4. Scattering from a zz{5757) and a 2:2^(558) defect line 

In this section we consider two other types of extended defect lines, namely the zz{5757) 
[16] and the zz{558) defect lines [TTl [H [7] [see Fig. [Tj. These defect lines are more 
realistic, albeit somewhat more complex to treat. We will see that both these defect 
lines give rise to a duplication of the unit cell. This originates a folding of the Brillouin 
zone, that brings to play two additional scattering modes. Despite the fact that at 
low-energies the latter happen to be evanescent modes, they must be considered in the 
computation of the scattering amplitudes. 



a 



(b) 



2ui 



2 \Mc 


k 
















1 


! 










\, jBi(n) TB^in) 



n = -1 



n 



m — 1 



m m + 1 




m + 1 



Figure 7. (Color online) Panel (a): Scheme of a zz(5757) defect line. Panel (b): 
Scheme of a zz(558) defect line. 



The treatments of both the zz{b7b7) and the zz{bb'8) defect lines are very similar to 
each other. The only difference between them is on the microscopic details of each defect 
line. As a consequence, and for the sake of definiteness, we shall first study electron 
scattering from a zz{b7b7) defect line. We will employ again the procedure described 
in Section [2| which was previously applied to the study of the pentagon-only defect line 
(see Section |3]). Afterwards, we will also compute the boundary condition arising from 
the 2;2;(558) defect line and we present the results for both the zz{b7'b7) and the zz{bb'8) 
cases. 



Formulating the problem 

Let us thus concentrate on the study of a graphene layer with a zz{b7b7) defect line 
along the zigzag direction. As can be seen in Fig. [7]^a), the zz{h7h7) defect line has a 
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periodicity twice as large as that of the pentagon-only defect hne [see Fig. |3j. Therefore, 
its unit cell will necessarily be two times bigger in the defect direction, when compared 
with the unit cell of pristine graphene. 

Then, the most natural choice for the lattice vectors is 2ui and U2 [where ui and U2 



were defined in Eq. (15)]. As a consequence the unit cell of bulk graphene will be twice 
as big in the Ui direction. Therefore, the First Brillouin Zone (FBZ) will be folded in 
comparison with the one of pristine graphene (see Fig. |8]). The Dirac points will then 
be located at Kj, = z/7r/(3a)(l, —y/S) (see Appendix D), where again, u = ±1 identifies 
the Dirac point. 



(a) 








y^,^ 




\ ■ 









Figure 8. (Color online) Panel (a): Scheme of the FBZ of pristine graphene. Panel 
(b): Scheme of the FBZ of graphene with a doubled unit cell in Ui direction. 

Let us now write the TB Hamiltonian of the graphene layer with a zz{5757) defect 
line [see Fig. [7]^a)]. As before, we can separate the Hamiltonian in three distinct parts: 
H = + + H^, where {H^) stands for the Hamiltonian above (below) the 
defect line, while the remaining term, , describes the defect line itself. In second 
quantization the explicit forms of and read 

ffU{L) ^ -t^^ (b\{m,n) + b\{m,n-l) + bl{m-l,n-l)^ai{m,n) 



+ (^bl{m, n) + foj (m, n — 1) + b\{m, n — 1) ja2(m, n) + h.c. 



(34) 



where for (H^) for n > 1 (n < —1). 

Similarly, the term describing the TB Hamiltonian at the defect line reads 



H 



D 



- (^^cd\{m) + b\{m, -1) + bl{m - 1, -l)jci(m) 

+ (^c4M + b\{m, -1) + bl{m, -l))c2(m) + ^ad\{m)d2{m - 1) 
+ UKm)f2{m) + ^,d\{m)f,{m) + ^J,{m)f2{m) + (^UKm) 
+ a\{m, 1) + 4(m, l))^i(m) + (^UIM + a\{m + 1, 1) 
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+ al(m,l)^g2{m) + h.c. , (35) 

where C,a, and C,c stand for the renormahzations of the hopping amphtudes at the 
defect [see Fig. [7](a)]. 

Again the system is invariant under translations r = 2mui (where m is an integer). 
Consequently, we can make use of Bloch theorem and Fourier transform the Hamiltonian 
along this direction, diagonalizing it with respect to the variable m. As expected, the 
resulting Hamiltonian turns out to be equivalent to the Hamiltonian describing a quasi- 
one-dimensional chain, with four atoms per unit cell and a defect at its center (see Fig. 
9l). 




n = -2 



n 



n = 



n = 1 



n = 2 



U2 



Figure 9. (Color online) Effective one-dimensional chain obtained after Fourier 
transforming the Hamiltonian of graphene with a zz(5757) defect line. The hopping 
parameters t' and ?/, are short hands for t' = te~^'^^''°' and for -q = t£^ae~^'^^''°' ■ The 
parameters ^a, Cfc ^-^d stand for the hopping renormahzations appearing in Fig. [vj^a). 

Therefore, we can write the Hamiltonian of the effective one-dimensional chain as 

H{kx) = H^{kx) + H^{kx) + H^{kx). The part above and below the defect read 

H''^'^\k,)= -tJ2 {b\{k,,n) + b\{k,,n-l) + e'^'^%l{k,,n-l)^ 

n '- 

X ai{k^,n) + {b\{k^,n) + b\{k^,n- 1) + b\{k^,n- 1)) 



X a2{kx, n) + h.c. 
while the part corresponding to the defect reads 



(36) 



icd\{k^) + b\{k^, -1) + e'^''^''b\{k^, -l))ci{k,) + (icd\{k^) 

+ b\{k^,-l) + b\{k^,-l)y2{k^) + U-'^^^'^d\{k,)d2{k,) 
+ Ul{K)f2{k,) + iA{k^)h%) + iA{k.)f2{k^) + [uKK) 
+ a\{k,, 1) + 4(A:,, l)]gi{h^ + [uliK) + e'^^^^'^aUA;,, 1) 
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+ al{h, l))g2{h) + h.c. 



(37) 



As was stated in Section [STTI the TB equations of this system are obtained from its 
Schrodinger equation, H{kx)\fi, k^) = e^^kJfJ',kx)- In this last equation, H{kx) stands 
for the Hamiltonian of the effective one-dimensional chain given by Eqs. (p6|-([37|). At 
the same time, the eigenstate k^) can be expressed as a linear combination of the 
site amplitudes along the one-dimensional chain: 

2 

\^i,kx) = X^X] |^Ai(fc^,n)|ai; fc^,n) + Bi{kx,n)\hi] k^.n) 

i=l n^O 
2 



+ ^Ci{kx)\ci; kx) + Di{kx)\di] k^) + Fi{kx)\fi] k^) 
1=1 



+ Gi{kx)\gi] kx) 



(38) 



where the \zi] kx, n) = zj^kx, n)\0) stand for the one-particle states at the atom Zi of unit 
cell n of the one- dimensional chain [with z = a,b on the bulk (n ^ 0) and z = c,d, f,g 
at the defect (n = 0)]. 



4-2. Bulk properties 

From the above discussion, we can readily write the TB equations in the bulk (n ^ 0) 
of the quasi-one-dimensional chain. They read 

eAi{kx,n)= - tBi{kx,n) - tBi{kx,n - 1) - te-'^''''''B2{kx,n ~ 1) , (39a) 

eBi{kx,n)= -tAi{kx,n)-tAi{kx,n + l)-tA2ikx,n + l), (39b) 

eA2{kx,n) = - tB2{kx,n) - tBi{kx,n - 1) - tB2{kx,n - 1) , (39c) 

eB2{kx,n) = - tA2{kx,n) - te'^''^''Ai{kx,n + 1) - tA2{kx,n + 1) . (39d) 

As before, we can write a recurrence relation between the amplitudes of the unit 
cell located at position n and those of the unit cell located at position n + 1. Following 



a procedure similar to that leading to Eqs. (23) and (24), one obtains 

L{n + l) = T{e,kx)Mn), (40) 

where L(?7,) = [Ai{kx,n), Bi{kx,n), A2{kx,n), B2{kx,n)]'^ , and the transfer matrix, 
T(e, kx), has the following form 



T(e, kx) 



oi2kxa 



-1 



1 



— e 



-e w{e,kx) 

-w{e-kx) 



-w{e, -kx 
w(e, kx) 



(41) 



— \-\- g*2fci:a _j_ g2 



where w{e,kx 

If we take into account the nature of the folding of the FBZ, we can make a basis 
change that will turn the transfer matrix T(e, kx) into a block diagonal matrix (see 



Appendix C). Such a basis change groups the modes in two distinct pairs, uncoupling 
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the bulk descriptions of each one of these pairs. We dub the modes of one of those pairs 
as the + modes, while the modes of the other are going to be called — modes. Each 
one of these pairs of modes is associated with a pair of energy bands present in the 
folded FBZ. Note, for instance, that when we position ourselves near the Dirac points, 
Kj, = i/7r/3(l/a, — a/S/o), the + modes turn out to be high in energy (e+ ~ 2t), while 
the — modes are low in energy (e_ ~ 0) |18j . 
The latter basis change reads 

Lin)=A{K)-L{n), (42) 

where L(n) = [A^{kj.,n), B^{kx,n), A_{kx,n), B^{kx,n)]^ , while the matrix A{kx), 



reads (see Appendix C) 



1 

71 





1 



1 



-ikxa 











(43) 



Applying this transformation to the transfer matrix (41) we obtain T{e^kx) 
A(A;,)T(e, A;,)A-i(A;,), where 



T(e, K, 



(44) 



T+(e,A;.) 
T_(e,A;,) 

is block diagonal. In the above expression, and T_ are 2x2 transfer matrices 
associated with the + and the — modes, respectively. These matrices are given by: 

1 



' 2 



2 cos 



> 2 ^ 



-f 4 COS 



2 ^ fc^a^ 



t2 



kxa) 



' 2 



2i sin 



- 2 > 



4sin^( 



t 



t2 



(45a) 



(45b) 



Again, as was said in Section [2| the eigenvalues and eigenvectors of the transfer 
matrix T give the Bloch solutions of the problem. Eigenvalues with unit modulus, 
|Ap = 1, correspond to Bloch solutions propagating along the one-dimensional chain (a 
band state), while eigenvalues with non-unit modulus, |Ap 7^ 1, correspond to evanescent 
Bloch solutions [decreasing with n — )■ -|-oo {n — )■ —00) when |Ap < 1 (|Ap > 1)]. 

Moreover, in the basis that uncouples -|- and — modes the matrices and T_ 
give the Bloch solutions associated with each one of the -|- and — energy sectors. Again, 
these can be either propagating or evanescent depending on the energy e and momentum 
kx- Of relevance is the case where the scattering process occurs at a sufficiently small 
energy (e < 2t) so that we are sufficiently near the Dirac points. In this case, the -|- 
modes turn out to be evanescent while the — ones are still propagating. When we are 
at very low-energies, e ~ 0, the transfer matrices associated with each pair of modes 
read 





T+ 0,z/ 



TT 



3a 



— e 



1 

V3 



(46a) 
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T_ 0,z/- 



TT 

3a 



1 
1 



(46b) 



Thus, we can conclude that the + modes are evanescent either decreasing or increasing 
exponentially as e~"^°s^ or as e"^°§^ with increasing n. At the same time the — modes 
are propagating. 



4.3. The zz (5757) defect 



It is important to note that the conclusions of Section 4^ apply not only to the case 
of the zz{5757) defect line, but also to the 2:2(558) defect line. In fact, the reasoning 



pursued in Section |42] describes the bulk of systems whose unit cell's size in the direction 
of Ui is twice that of pristine graphene. As we can readily conclude from Fig. [7], both 
the 2;2;(5757) defect line and the zz{558) one impose the same duplication of the unit 
cell along the zigzag direction. As a consequence, the bulk problem of these two systems 



is described by the same transfer matrix relations, Eqs. (41)-(45). 

As previously noted, what distinguishes electron scattering in the zz{5757) and in 
the zz{558) defect lines are the microscopic details of each one of these defects. Next 
we will show how to compute the boundary condition relation for each one of these two 
defect lines. 



Let us start by the case of the 2;2;(5757) defect. From Eqs. (37) and (38) we write 
the TB equations at the 2;2;(5757) defect. They read 



- ^A(l) 


= W\G{0) + B{1), 


(47a) 


^G(O) 


= ^,F{0) + WaA{1), 


(47b) 


^F(O) 


= e6D(0)+ecG(0)+eafx.F(0), 


(47c) 


^D(O) 


= ecC(0)+e6F(0)+ea<D(0), 


(47d) 


^C(O) 




(47e) 


^B(-l) 


= A{-1) + WaC{0), 


(47f) 



where Z(ra) = [Zi{n), Z2{n)]'^ , for Z = A, B,C, D, F,G. The matrices Wa, cTx and a'^ 



appearing in Eqs. (47|), are given by 
Wa = 



1 



1 

1 



e 



(48a) 



(48b) 



(48c) 



In close analogy with what was done for the case of the pentagon- only defect line, 
and after some elementary algebraic manipulations, we can rewrite the above 2x2 
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matrix equations in a more compact form: 



B(l) 
A(l) 

A(l) " 
G(0)_ 

G(0)" 
F(0) _ 

F(0) " 
D(0) _ 

D(0)" 
C(0) _ 

C(0) 
B(-l) 



fl2 







-I2 



A(l) 
G(0) 



-I2 



4 




G(0) 
F(0) _ 

" F(0) 
_D(0) 

"g(o) 

C(0) 



-h 
HWa)- 







C(0) 
Bf-11 



(Wa) 




B(-l) 
A(-l) 



(49a) 
(49b) 
(49c) 
(49d) 
(49e) 
(49f) 



where I2 stands for the 2x2 unit matrix. Note that Eqs. (49) are now 4x4 matrix 
equations. In what follows we are going to denote the matrices in Eqs. (49a), (49b), 
( F9^ , ( |49dl ), (ig and (|4§) by Pi(e,fc,), F^ie,^), Fsie,^), F^{e,k,), F^ie,^) and 
^6(^5 kx), respectively. 

The boundary condition matrix equation relating the vectors 

L(l) = [Aiih, 1), B,ih, 1), A2{h, 1), B2ih, l)f 



and 



becomes 



L(-l) = [A(fc,, -1), B,{K, -l),A2iK, -1), B^iK, 



(50) 
(51) 



(52) 

5757, is a 4 X 4 matrix (in 



L(l) =M5757.L(-1). 

In the above equation, the boundary condition matrix 
contrast with the case of the pentagon-only defect line) and is given by 

M5757 = i?.Pl.P2.P3.F4.P5.P6.i?^, (53) 

where, for the sake of simplicity of notation, we have omitted the dependence of the 
matrices Pj on e and on k^; the matrix R in Eq. (52) makes a basis change from 
{Bi, B2, Ai, A2} to {Ai, Bi, A2, B2}, and is defined by 

10 



R 



(54) 



10 
1 
10 

It is worth commenting that, if we change to the basis that uncouples the + and 
energy sectors, M5757(e, fc^;) = AM5757(e, /c2^)A~-'^ [where matrix A is defined in Eq. 



(43)], we readily conclude that the boundary condition matrix mixes + and — modes 
on opposite sides of the defect. 
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We can employ an entirely analogous procedure to the zz{558) defect. As before, we can 
write the Hamiltonian describing such a system, and from it we write the TB equations 
at the zz{558) defect [see Fig. [7]^b)]. These read 



-^B(O) 
-^D(O) 

-!a(o) 



t 



B( 



W\B{0) + B{1), 
aD(0) + W^^A(l), 
6B(0)+6<D(0)+6A(0), 
6D(0) + W^iB(-l), 
A{-1) + WaA{0), 



(55a) 
(55b) 
(55c) 
(55d) 
(55e) 



where, once more we use the notation 7i{n) = [Zi{n), Z2{n)Y , for Z = A,B, D. 

As before, we can easily rewrite the above 2x2 matrix equations in a more compact 

form 



B(l) 
A(l) 

A(l) " 
B(0) _ 

B(0) " 
D(0) _ 

D(0)" 
A(0) _ 

A(0) 
B(-l) 







A(l) 
B(0) 



-h 



-I2 

iiWA) 





I2 





B(0) 
D(0) 

D(0) 
A(0) 



-h 

HWa) 
-h 



A(0) 

Bf-r 



-1 



{Wa 




,-1 



B(-l) 
Af-11 



(56a) 
(56b) 
(56c) 
(56d) 
(56e) 



The 4x4 matrices appearing on the right hand side of Eqs. (56a), (56b), (56c), (56d) 
and ( |56e| ) will be denoted by Qi{e,k^), Q2{e,K), Q3{(^,kx), Qi{(^,kx) and Q5(e, /^x), 
respectively. 

The boundary condition matrix equation relating L(l) and L(— 1) is of the same 
form as in Eq. ( [52) ) 

L(1)=M558.L(-1), (57) 

where, the boundary condition matrix, M558, is a 4 x 4 matrix now given by 

M558 = R-Qi-Qi-Qs-QA-Q^.R" , (58) 

where again, for the sake of simplicity of notation, we have omitted the dependence of 
the matrices Qj on e and on k^- 

Once again, in the basis that uncouples the + and the — modes, M558(e, A;^) = 
AMssgA"^ [where matrix A is defined in Eq. (43)], one readily sees that the boundary 
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condition matrix also mixes the modes of the + and the — energy sectors of opposite 
sides of the defect. 



4-5. Computing the scattering coefficients 



The scattering problem from the zz{5757) [and the zz{558)] defect line can now be 



solved completely using the boundary condition at the defect, Eqs. (53) [Eq. (58)], and 



the modes allowed in the bulk [obtained from the transfer matrix, Eq. ([44j)-([45|)]. 

Let us suppose that we are working at positive energy and around the u = +1 
Dirac point [K+ = 7r/3(l/a, —Vs/a)]. Moreover, we suppose that the energy is small, 
so that the + modes are evanescent. We then consider the scattering process in which 
we have one incoming mode from n = — oo associated with the — energy sector. Due 
to the presence of the defect at n = 0, there will be a reflected as well as a transmitted 
mode. Besides, in accordance with what was said above, there will be two evanescent 
modes associated with the + energy sector. 

We choose to denote I^E'^^) (|\1'1^^)) as the right (left) moving mode associated 
with the — energy sector of the transfer matrix [see Eqs. (41)-(45)] when we are both 
at e > and around the K+ Dirac point; the corresponding eigenvalues are denoted 



It') 



Similarly, we denote by |\&^'*) (|\E'^'')) the + energy sector's transfer 



(+)\ 



by A^-) (A 

matrix mode [see Eq. (41)-(45)], which decreases (increases) in the direction of U2. This 
mode's corresponding eigenvalue is going to be denoted by A^^^ {X[^K) We can thus 
write the wave-function on each side of the defect as 



L(n < 0) = (a1:^)"+'|^^")) +p_(Al7))"+'|^lr^) 



L{n > 0) 



(59a) 
(59b) 



where p_ and r_ (p+ and r+) are the reflection and transmission scattering amplitudes 
of the — (+) modes, respectively. 

As before, for different choices of e and kx the incoming, reflected, transmitted and 
evanescent modes must be chosen accordingly with the direction of propagation and/or 



direction of increase of the eigenmodes of the transfer matrix (see Fig. B2 ) 



In order to determine the scattering coefficients p± and t±, and in accordance with 
Section |2| we just need to substitute Eqs. (59) in the boundary condition expression 
arising from the zz{5757) [zz{558)] defect, Eq. (52) [Eq. (57)]. to obtain the equivalent 



of Eq. ( 12 ) for these cases. 





r_ 




U-^MU 



J n=+l 



P- 



(60) 



where we have omitted the subscripts identifying the boundary condition matrix [either 
M5757 for the zz{5757) defect, or M558 for the zz{558)]. In the above equation, the matrix 
U is the matrix mediating the basis transformation from the basis {A^, B^, A_, B_} to 
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the proper basis of the transfer matrix. As stated in Section [2| each column of U is one 



vi/W),|vi/^-o,i^r),i^ 



,(+) 



(-)n 



i:'} 



written in 



of the eigenmodes of the transfer matrix U = 
the 5+, basis. 

Solving the above linear system of four equations can be readily accomplished with a 
computer algebra system. Its solution determines completely the scattering amplitudes 
of the problem. 



4.6. The transmittance 

As before, we are interested in computing the transmittance for relatively low-energies, 
which is equivalent to say that we want to compute the transmittance around the 
Dirac points. Consequently, the + modes will be evanescent while the — modes will be 
propagating. Therefore the transmittance is going to be given by T = jr. p. 

Similarly to what was done for the transmittance in the pentagon- only defect case, 
here, and for the sake of comparison between scattering processes occurring at positive 
and negative energy, we will plot the transmittance against 6' = 6 when e > 0, plotting 
it against 6' = —6 when e < 0. 



In Fig. 10 we represent the transmittance T = |r_p [associated both with the 
^^(5757) and the zz{558) defect lines] as function of the angle 9 (when e > and —9 
when e < 0), made between q and the defect line. In the several panels of Fig. 10, we 
present the transmittance curves for different values of the energy [the ones referring to 
the zz{5757) defect are those in panels (a), while the ones referring to the zz{558) defect 



are those in panels (b)]. The curves in the Fig. 10 refer to the Dirac point K4.. Those 
referring to the Dirac point K_, are mirror symmetric relatively to the line = 7r/2 to 
the ones presented in the figure. In order to obtain the transmittance curves of Fig. [10} 
all hopping renormalizations were set equal to one in both defects: = ^fe = = 1 and 

6 = 6 = 1. 

As before, modifications of the hopping parameters at either the zz{5757) defect 
of the 2:2; (558) defect, result in very different transmittance behaviors. Moreover, the 
low-energy result is in good accordance with that obtained in the continuum limit in 
Ref. 



5. Conclusions 



In this paper we have presented a procedure to work out the electronic scattering either 
from a defect in a one- dimensional crystal, or from a periodic extended defect in two- 
dimensional crystals, which can be reduced to a quasi-one-dimensional tight-binding 
model. The formalism developed only uses simple tight-binding concepts, which reduce 
the scattering problem to a set of matrix manipulations. These can be easily worked 
out by any computational algebraic calculator. 

We have illustrated the presented procedure in the context of the first neighbor 
tight-binding model of graphene when a defect line is oriented along the zigzag direction 
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Figure 10. (Color online) Panels (a): Transmittance through the zz(5757) defect 
line as function of the angle 6 (when e > and —6 when e < 0) between q and the 
barrier. The hopping parameters at the defect have the value £,a = £.b = £.c = ^- Panels 
(b): Transmittance through the zz(558) defect line as function of the angle 6 (when 
e > and —9 when e < 0) between q and the barrier. The hopping parameters at 
the defect are = ^2 = 1- In each group of panels, the transmittance is plotted for 
several different energy modulus. These are: (al) and (bl) |e/t| — 0.01; (a2) and (b2) 
|e/t| = 0.05; (a3) and (b3) \e/t\ = 0.1; (a4) and (b4) \e/t\ = 0.3. Positive energies 
are represented by the black full lines, while the negative ones are represented by the 
dashed red lines. Only the curves associated with the Dirac point K_|_ are represented. 
Those for the Dirac point K_ are obtained from the former by a reflection of these 
over the axis 9 — 7r/2. 

in an otherwise perfect crystal. Three distinct defects were studied: the pentagon-only, 
the zz{5757), and the zz{558) defect hnes. 

The latter two defect lines forced a duplication of the unit cell along the defect 
direction (relatively to the case of pentagon-only defect line). Such duplication renders 
the mathematical treatment of the problem more complex, by introducing an additional 
pair of scattering modes. These latter modes happen to be propagating at high 
energy around the Dirac points. Thus, for low-energy scattering processes they will 
be evanescent states. 

All the three defect lines studied have different behaviors that strongly depend on 
the particular values of the hopping parameters at the defect lines. Moreover, in a 
companion paper we treat these same problems in the continuum low-energy limit 
using the Dirac equation. We have shown that we can recover the low energy results 
from the tight-binding approach, as it should be the case. 

Noteworthy is the fact that the present procedure can be used to treat more realistic 
extended periodic defects, that do not need to be linear. In principle, such kind of 
extended defects will force a considerable increase in the size of the unit cell, thus 
introducing a larger number of scattering modes. This, of course, renders the problem 
a more difficult one. 
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Appendix A. Left and right eigenvectors of the transfer matrix T 



The right eigenvectors of a non-hermitian matrix T, are defined by the equation 

T\iJ^) = X^\^^), (A.l) 

while the left ones, are defined by 

{i^^\T = \Si\. (A.2) 
It is straightforward to demonstrate that, apart from a normalization constant, 



Ci = y {ipilipi), they form a dual basis. Explicitly, 

(V^,|T|^,) = Xj{ip,m = A,(V^,#,), (A.3) 

and then, if Aj 7^ Xj, we must have = 0. If we define \ipi) = \4'i)/Ci and 

\{pi) = \ipi)/Ci, we can thus write 

{^j\LPi) = Sij. (A.4) 

Let us now show that the left eigenvectors of a non-hermitian matrix, T, are the 



right eigenvectors of its transpose, T^. If we take the hermitian conjugate of Eq. (A.2), 
we obtain 

Tt|V^,) = A*|V^,). (A.5) 
Spelling out this equation in components, and taking the complex conjugate, we obtain 

= A.dv^.)):, (A.6) 

which gives us the components of the row vector (V^j|a = {\4'i)Ta- 
Appendix B. The conserved current along the ID chain 

In this appendix we are going to compute the conserved current along a quasi-one- 
dimensional chain. We start by doing it for a general case (using the formalism of 
Section [2]). Later, we write the conserved current for the one-dimensional chains arising 
from the Fourier transformation of both the graphene layer with a pentagon-only defect 
line, and that with a zz{5757) [or a zz{55S)] defect line. 

The bulk TB equations of the general quasi-one- dimensional chain, Eq. ([T]), can be 
rewritten in the following manner 

ec{n) = Vlc{n - 1) + Hrc{n) + VLc{n + 1). (B.l) 

where we have used Vr = Vl (Hamiltonian hermiticity) . The time-dependent 
counterpart of this equation can be written as 
d 

z^— c(n, t) = Vlc{n - l,t) + HrC{n, t) + VLc{n + l,t). (B.2) 



The adjoint of Eq. (B.2) reads 



d 

-ih—c^n,t) = c^{n-l,t)VL + c^{n,t)Hr + c^{n+l,t)Vl (B.3) 
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where 
write 

where 



Hr from the TB Hamiltonian hermiticity property. Given this, we can 
d 



dt 



\n,t)c{n,t)= - J{n,t) - J{n-l,t) 



(B.4) 



J{n,t) 



(B.5) 



: c^{n,t)VLc{n + l,t) -c^{n + l,t)Vlc{n,t) , 
ill J 

can be interpreted as the particle current along the one-dimensional chain flowing from 
position n to position n + 1. If we have a stationary state of an AB-\ike chain, we can 
use the transfer matrix relation [see Eq. (pi)] and write the conserved current as 



Jin) 



cUn) 



VrJ - TV/ 



cm 



(B.6) 



We can further rewrite the conserved current expression using the fact that Bloch modes 
are the eigenvectors of the transfer matrix, Tipx = Xipx. Therefore, for a Bloch solution 
labeled by A, whose eigenvector is ipx, the conserved current reads 



Jx 



XVr 



\*vl 



For the case of the pentagon- only defect line it can be inferred from Eqs. 
the matrix Vl reads 



1 + e^^-" 



Vr 



- 1 



(B.7) 
21 p that 

(B.8) 



while the transfer matrix is given by Eq. (24). As a consequence, we can easily conclude 
that the conserved current operator along the one-dimensional chain (originated from 
the Fourier transformation of pristine graphene along the direction of Ui), is given by 

' (B.9) 



J ='-{yLT-T^vl 



(B.IO) 



where ay is the y-Pauli matrix. As a consequence, the conserved current reads 
Jx = ^^{(^yi^x- 

Note that in the context of the two-dimensional graphene lattice, the above 
conserved current is a current directed along U2. Therefore, the current along the 
y-direction is given by Jy{X) = J'xy/3/2, which is in accordance with the low-energy 
result. 



In Fig. Bl we present a density plot of the current associated with each one of the 
modes of the transfer matrix of pristine graphene. 

Similarly, for a one-dimensional chain obtained from the Fourier transform of 
graphene with a doubled unit cell in the direction ui [which happens in the cases of 
the zz{5757) and zz{558) defect lines], we can infer from the bulk equations, Eqs. (39), 
that the matrix Vl reads 





Vr 



t 





1 



^2 Af ^ Q> 










(B.ll) 
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Figure Bl. (Color online) Panels (a) and (b) are density plots of the current 
associated with the modes of the transfer matrix of the pristine graphene lattice. In 



accordance with the definition of sub-Section 3.4 panel (a) gives the current associated 



with the mode |\E'>), while panel (b) gives the current associated with the mode |^'< 



while the transfer matrix for this quasi-one- dimensional chain is given by Eq. (41). As 



a consequence, we can readily conclude that the conserved current operator along this 
chain, is given by 



J 



h 



TV2 



t 

h 



Oy 

Oo 



09 



cr,; 



(B.12) 



where Oy stands for the |/-Pauli matrix, while O2 represents a 2 x 2 zero matrix. Note that 
the appearance of two copies of pristine graphene conserved current was to be expected 
because, as was already referred, the folding of the FBZ due to the duplication of the 
unit cell along the direction of ui brings an additional pair of modes into the folded 
Brillouin zone. Both the modes associated with the + and the — energy sector were 



already known to have a current operator given by Eq. (B.9). 



In Fig. B2 we present a density plot of the particle current associated with each 
one of the modes of the transfer matrix of graphene with a doubled unit cell along the 
direction of Ui. 



Appendix C. The basis uncoupling the high and low-energy modes 

As we have argued in the main text it is, in some circumstances, required to work with 
a doubled unit cell. That is, for instance, the case when treating the Z2;(5757) and the 
zzi^'b'S) defect lines in graphene. 

With either of these two defect lines present, Ui is no longer a lattice translation. 
To accommodate this, it is convenient to describe pristine graphene with a unit 
lattice defined by the vectors 2ui and U2, instead of the usual choice, ui and U2. 
The corresponding reciprocal lattice vectors will then be vi/2 and V2, instead of 
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(c) 




1.0 -0.5 0.0 0.5 1.0 1.5 




(d) 



e 





l.B -1.0 -0.5 0.0 0.5 1.0 1.5 



Figure B2. (Color online) Panels (a)-(d) are density plots of the current associated 
with the modes of the transfer matrix of the graphene lattice doubled along the Ui 



direction. In accordance with the definition of sub-Section 4.5 panel (a) [(b)] gives the 
current associated with the (+) modes, while panel (c) [(d)] gives the current associated 
with the (— ) modes. 



those associated to the FBZ of pristine graphene, vi = 27r/a(l, l/-\/3) and V2 = 
27r/a(0, 2/\/3). Therefore, the doubhng of the unit cell will have as its main consequence 
the folding of the FBZ along the direction of vi [compare both panels of Fig. [s]. 

There will thus be twice as many atoms in the doubled unit cell as those contained 
in the pristine graphene one [two atoms of sub-lattice A [Ai, A2) and two atoms of 
sub-lattice B {Bi, B2)]. At the same time, there will be twice as many energy bands 
in the folded FBZ as those present in the pristine graphene one. Moreover, the Dirac 
points will now be located at K.^ = z/7r/3(l/a, —^/3/a) (see Fig. |8|. Near the new Dirac 
points, two of the four energy bands present in the folded FBZ are the two low-energy 
Dirac cones, while the other two bands are the high-energy bands flUi [T5j . In what 
follows, we identify the two bands of low-energy near the Dirac points by — , while the 
high-energy ones are identified by -|-. 

The bulk physics of graphene with a doubled unit cell must be exactly the same 
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as that of pristine graphene. In fact, the Bloch vectors q and q + vi/2 identifying 
two different Bloch wave solutions in the unfolded Brillouin zone of pristine graphene 
correspond to a single Bloch vector q in the folded zone. As a consequence, in the case 
of graphene with a doubled unit cell it is natural to expect that we can find a basis in 
which we can uncouple the physics associated with each one of the two Bloch solutions 
of the unfolded system: that identified by the Bloch vector q and that identified by the 
Bloch vector q + Vi/2. 

Let us then put ourselves at the wave vector k = q, where q is close to the Dirac 
point of the folded Brillouin zone. Bloch theorem in the unfolded pristine graphene 
allows us to write the following relations between amplitudes of the doubled unit cell: 

^2 = e^'i-^Mi = e'^-Mi, (C.la) 
B2 = e'^'-'^'Bi = e'^^^Bi. (C.lb) 

Similarly, for the wave vector k = q + vi/2, we have the relations 

A2 = e*('i+^)-"Mi = -e'^^^'Ai, (C.2a) 
^2 = e*('i+^)-"^fii = -e*'?-''5i. (C.2b) 

We note that the vector q + vi/2 can be close to the Dirac point of the unfolded zone. 
In fact, if we choose q = = (vr/Sa, —n/\/3a), the Dirac point of the folded zone, 
q + vi/2 = (47r/3a,0), the Dirac point of the unfolded zone. Then, in the unfolded 
zone, the state with k = q corresponds to a high energy state (+), whereas the state 
k = q + vi/2 corresponds to a low energy state (-). This analysis motivates the past 
definitions of high and low energy modes. 

As said above, we want to construct a basis that uncouples the modes originating 
at the two different locations of pristine graphene's FBZ: k = q and k = q + Vi/2. 
Equivalently, we want to construct a basis that verifies A^ 7^ 0, 5+ 7^ and 
A_ = = when k = q, as well as, A^^ = B^ = 0, A_ 7^ and B^ 7^ when 
k = q + vi/2. Given this, we define the new basis as 

^+ = ^(A + e-^''^M2), (C.3a) 
v2 

fi+ = -^(Ai + e-'^^'^52), (C.3b) 

A„ = ^(Ai-e-^'^-Ms), (C.3c) 
v2 

B_ = ^(Ai - e-*^-"52). (C.3d) 
V2 

which verifies the previous conditions. This new basis uncouples the (+) and the (— ) 

energy sectors of the doubled unit cell. Note that the new basis written in Eqs. (C.3) 
is exactly the same defined in Eqs. (42) and (43). 
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Appendix D. The Dirac angle 9 



In the main text the transmittance is given in terms of e and 0. On the other hand, in 
the figures appearing in the main text the transmittance is given as function of e and 9 



[see. Fig. Dl for the definition of 9]. We show now how given e and we can obtain 9. 
We consider the case of a doubled unit cell. 

For a doubled unit cell Schrdinger's equation reads 



E 



' A, ' 
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h 





t2 ' 






Bi 




f* 












Bi 


A2 












h 




A2 


B2 




''2 





''I 







B2 



(D.l) 



where 



t2 



-til + e 



-ik u2 ^ 



(D.2) 
(D.3) 



Transforming the Hamiltonian to the block diagonal basis we obtain 



E 



A+ ' 







i+ 










' A^' 


B+ 




t* 













B+ 


A_ 













t_ 




A_ 


B_ 

















5_ 



(D.4) 



where 



til 



'i{k-U2+<l>)) 



- t(l + e-'^-'"' - e-'^''-''^- 



(D.5) 
(D.6) 



We can find a set of momentum values where t_ is equal to zero. This corresponds to 
the new position of the Dirac points in the new Brillouin zone. The quantity t- can be 
written as 



t- = -t(l + e 



-i(V3kya/2~4,/2) _ ~i{^kya/2+4>/2) 



) 



t(^l + 2ie-''^^y''/^ sin 



If sin(A;x.a/2) = ±1/2 and e 

TT 



-i\/3kya/2 



±i we have t_ = (and t^ 



k^a = ±- 



JvyCL 



TC 



Thus, the new Dirac points are located at 
Kd-- 



^TT TT 

3a' a/So 



(D.7) 

-2t). This implies 
(D.8) 

(D.9) 
(D.IO) 



The Dirac angle 9 is defined as in Fig. Dl 



The dispersion of the low energy modes reads 



\t-\Vt' 



1 + 4 sin — 1-4 sin — sin 



ky\/3a 



(D.ll) 
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Figure Dl. (color online) Definition of the Dirac angle 9 € [0, tt]. 



The largest and the smallest value of = ak^ are given when ky = K£,y. In this case 



2 sin 



lie 



sm 



Then 



and 



2 arcsin 



1 -e 



0_ = 2arcsm — 
The coordinate kya is given by 

/c,,a = arcsin — 
^ 73 



4sin^ 



4 sin 



Since k = K^, + q we have 

Qx — kx -fi^Dx = 



Qy — ky -^Dy 



ky ~\~ 



From Fig. Dl is clear that 



9 = arccos 



TT 

3a ' 
a\/3 



arccos 



0-7r/3 



(D.12) 
(D.13) 
(D.14) 

(D.15) 

(D.16) 
(D.17) 

(D.18) 



Thus, given and e we can compute kya and d. Scanning between 
originates Q e [OjTt]. Thus ^ is a parametric function of and e. 



^min and (j)rnax 
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